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;> Abstract 

Q\ . We study the decay width and CP-asymmetry of the inclusive process b — ► sgg (g 

denotes gluon) in the three and two Higgs doublet models with complex Yukawa couplings. 
We analyse the dependencies of the differential decay width and C-P-asymmetry to the s- 
quark energy E s and CP violating parameter 9. We observe that there exist a considerable 
enhancement in the decay width and CP asymmetry is at the order of 1CP 2 . Further, it 



is possible to predict the sign of Cj using the CP asymmetry. 
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1 Introduction 

Rare B- meson decays, induced by flavor changing neutral currents (FCNC), involve at the 
loop level in the standard model (SM), therefore they are phenomenologically rich. The mea- 
surements of the branching ratio {Br), CP asymmetry (Aqp), forward-backward asymmetry, 
polarization effects, etc., provide restrictions on the SM parameters, such as the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants, etc. In addition, the 
possibility of replacing the SM particles by non-standard ones results in the sensitivity of these 
decays beyond the SM, like multi-Higgs doublet models, minimal supersymmetric extension of 
the SM (MSSM) 0, etc. The experimental effort at SLAC, KEK B-factories, HERA-B and 
possible future accelerators [0, [| stimulate the theoretical studies on these rare decays. 

Among B-meson decay modes, inclusive b —> sg decay becomes attractive since it is theo- 
retically clean and affected by loop contributions due to new physics beyond the SM. In the 
literature, there are various theoretical calculations on the Br of this decay. In the SM, the 
Branching ratio of b —>■ sg decay was calculated as Br(b — ► sg) ~ 0.2% for on-shell gluon 
||. However, to decrease the averaged cham multiplicity r\ c and to increase kaon yields || 
the enhancement of Brib — > sg) is helpful. The possibilities for the enhancement is the addi- 
tion of the QCD corrections and non-standard effects coming from the new physics. In |7|, [|, 
Br (b -► sg) was calculated in the 2HDM (Model I and II) for m H ± ~ 200 GeV and tan (3 ~ 5 
and it was found that there was an enhancement less than one order of magnitude. Further, 
this decay was studied in the supersymmetric models |§ and in the framework of model III 
2HDM [|Kj. In the model III, the enhancement was found at least one order of magnitude 
larger compared to the SM one and it was observed that there was no contradiction with the 
CLEO data [|TJ 

Br (b -> sg*) < 6.8 % (1) 

for light-like gluon case. Recently, 0(a s ) virtiual corrections and additional 0(a s ) brem- 



strahlung effects to the decay width of b — > sg was calculated in the SM |12[ and the en- 
hancement in the Br was obtained as more than a factor of two larger of the previous SM 
results. 

The inclusive process b — > sgg is another decay which has Br at the same order as Br (b — ► 
sg) according to the studies in the literature [T3|, T3, [T5|]. This process becomes not only from 



the chain decay b — > sg* followed by g* — > gg but also from the emission of on-shell gluons from 



the quark lines to obey gauge invariance. In [14], the complete calculation was done in the 



collinear and infrared singularity free region, in the SM and Br ratio was found at the order 
of magnitude 10~ 3 . In [l^, pjl the additional contribution of gluon penguins in the model III 
was estimated as negligible. Recently b — > sgg was studied in the model III with real Yukawa 



couplings [16] and a considerable enhancement was observed for the Br of the process even 2 
orders of magnitude larger compared to the SM case. 

In this work, we study the decay width T and the CP asymmetry Acp of b — > sgg decay 
in the model III and the 3if-DM(0 2 )- The reason to study the b —>■ sgg process is the possible 
considerable enhancement of T compared to the one in the SM and the measurable Aqp in 
the framework of the models underconsideration. In our theoretical calculations we choose the 



collinear and infrared singularity free kinematical region, following the procedure given in ][L4 
Here we take the source of CP violation as the complex Yukawa couplings appearing in both 
models. 

The paper is organized as follows: In Section 2, we give a brief summary of the 3HDM(02) 
and present the calculation of the decay width of the inclusive b — > sgg decay in the framework 
of the 3HDM(02) and the model III. Further we calculate the (differential) CP-asymmetry 
(Acp(E s )) A C p of the process. Section 3 is devoted to discussion and our conclusions. In 
Appendixes, we give the form factors appearing in the matrix element of the decay undercon- 
sideration and summarize the theoretical results for the 3HDM(02)- 

2 The inclusive process b — ► sgg in the framework of the 
model III and 3HDM{0 2 ) 

The general Yukawa interaction in 3HDM is 

£y = vfjQiLfaUjR + V%QiL(i>iD jR + tfjQiLfoUjR + gjQiLfaDjR 

+ pfjQiLfoUjR + pZQiLfoDjR + h.c. , (2) 

where L and R denote chiral projections L(R) = 1/2(1 =p 75), 0, for % = 1,2, 3, are three scalar 
doublets and 77^' , £y' , p t - are the Yukawa matrices having complex entries, in general. With 
the choice of scalar Higgs doublets 



1 

V2 



\ , ( V2V 



(3) 
, _ J_ / V2H+ \ _ J_ ( V2F+ \ 



and the vacuum expectation values, 

< 0i >= -7= f v J ; < 02 >= ; < 03 >= , (4) 

the SM particles carried by the first doublet and the information about the new physics by the 
others. The Yukawa interaction 

£y,fc = &jQiL<j>2U jR + ^QiLhDjR + pfjQiLfoUjR + pfjQiLfoDjR + h.c. . (5) 

describes the Flavor Changing (FC) one beyond the SM. Here, the couplings £, U,D and p u ' D for 
the charged FC interactions are 

£ch = 6v V C KM , 
ich — VCKM £n , 

Pch — Pn V C km , 

Pch = V C KM PN , (6) 

and 

C U,D _ fj/U,D,-l c U,D t/U,D 
SAT — \ V L ) S V R ■> 

v, D = {v u,o rY ,n v v, D (?) 

where the index "TV" in £ N ' denotes the word "neutral". Note that the Yukawa interaction 
for the model III can be obtained by taking into account only two doublets 0i, 2 and Yukawa 
couplings rjY j: £g. 

The decay amplitude of the process b — > sgg is given by 

M(b - sgg) = i^eZfaWkJsWTtfbb) , (8) 

where e%(k) are polarization vectors of the gluons with color a and momentum k, 

T ab _ T *_ ±_ , T E ±_ ±_ (q) 

± ixv ± p.v 2 o > IV 9 2 ' 

Here 4r are the Gell-Mann matrices and T/f, can be obtained by the replacements k\ <-> ki-, 
[A <-> v in the function T^. 

Since the process occurs at least at one-loop level in the SM, the function T^ v have con- 
tributions coming from light and heavy internal quarks, namely, u, c, t. Internal quarks d, s, b 
can also give contribution to the process beyond the SM. In the case of heavy internal quark, 



t-quark, the terms k 2 xternal /m 2 (rriyy, m H±, wipi) are neglected. However, for light internal 
quarks, k 2 xternal /m 2 terms can give considerable contribution. This forces us to parametrize 



the function T^ as 

T^ = T^ avy + T l ;f\ (10) 

(for the explicit forms of T^ avy and Tj^ ht see Appendix A). On the other hand, T" 6 can be 
divided into color symmetric and antisymmetric parts as JL4|] 

T$ = T+{±-} + IZ,[-,-], (11) 



with 



2 ' 2 J ^ L 2 



^ fJ,U ~ q I "'/•"' + Jyrjfc,, 



2 

^ = ^-T*). (12) 



Finally, using the expression 



T Sy m (Asym) _ Tr (C + (_) T + H (^ + m& )) p+H /) pW*' p™' ; ( 13 ) 

with the color factors C + = c 2jv c ~ anc ^ ^- = " 2 ~ an< ^ ^he polarization sum of the 
on-shell gluons 

pW _ — flW 4 _|_ _ X - ] 



we get the differential decay width of the process 

d 2 T 11 



\M\ 2 . (14) 



dE s dE t 27r 3 8m 6 ' 

Here E s is the s-quark energy and E\ is the energy of gluon with polarization e a {k\). M is the 
average decay amplitude, M = 2 j+\ w ^> anc ^ ^ = 2' ^ c = ^- 

In the expressions, the symmetric and antisymmetric parts do not mix each other. Further, 
the decay width can be divided into three sectors (see Appendix A): left (r L ), right {T R ) and 
left-right (T LR ). Left one is created by the nonvanishing k 2 xternal /m 2 ight terms, however right 
sector contains the forms factors with parameters mf/m^ and m 2 /m 2 H where % — u,c,t and H 
is one of the Higgs bosons. Left-right sector contains mixed terms. 

Now we are ready to calculate the CP-violating asymmetry Acp of the process b — > sgg. 
In the model III and 3HDM(02), the complex Yukawa couplings are possible sources for CP 



violation. Our procedure is to neglect neutral boson effects and all Yukawa couplings except 
gtftt and £$,bb (^N,tt ari d ^Nbb) ( see Discussion) in the model III (3S-DM(0 2 )) (see Appendix 
B for the definitions of e^ tt and e^ feb ). Therefore, in the model III (3HDM(0 2 )), only the 
combination ^Ntt^Nbb (^Ntt^Nbb) * s responsible for A C p. Using the parametrization 

1 &,«$,» e" (model III), 

4,« £nm (cos 2 + isin 2 6) (3HDM(Q 2 )) , 



m t m b 

1 7 V F D 



and the definition of differential CP asymmetry Acp(E s ) 



d 2 T fu > „„„\ d 2 T 



-(b -^ sgg) - m7m - i (b ^ sgg) 



a frp \ _ dE a dE-L *._ ; JJ ' ati 3 av\ ^ ^^" /,r\ 

dffES b - s ^) + d£fe(^ - S99)' 



_ d 2 r ^ , „„^ i d 2 v 

d 

we get 



AME.) = 2/ m (A„) G 2 fo) (*"(**■■- * *■))" (16) 



where 



fi = 18 Si m 6 (2 SI -m b )(-2S s + m b ) ((2 S s -m b )m 6 + 2 Si (2 S s + m b )) , 

A = -2 |S 2 | 2 m b (2496 ^ + 192 E\ (20 S s - 23 m b ) + mjj (-28S 2 + 44S s m b - 15 m 2 ) 

+ 2 E x m\ (172 S 2 - 208 E s m b + 69 m 2 ) - 4 S 2 m fe (316 S 2 - 562 S s m b + 213 m 2 ) 

+ 8 S x 3 (204 S s 2 - 638 S s 7^ + 357?/^)) 

+ 2 i?e(F 2 ) Se(A F 1 -Ai 2 ) c l + 8 E\ (2 S x - m b ) (-2 S s + m 6 )(7 | A i 5 \ 2 E s m b + 9 \Ai 2 \ 2 
(8 E\ + SE 1 (E s - m b ) + m fe (-3 S s + 2 m b )) 

+ 9 | A Si| 2 (8 S 2 + 8 Si (S s - m b ) + m fe (-3 E s + 2 m 6 )) 

- 18 J Re(AS 1 *Ai 2 )(8(S 2 H-SiS s -S 1 ra 6 ) -3S s m fe + 2m 2 ) . (17) 

Here 9 is the CP violating parameter which is restricted by the experimental upper limit of 
the neutron electric dipole moment eq. ( |22|) and F 2 = F^ HDM — S 2 (0), A Si, Ai 2 , and Ai 5 
are the Wilson coefficients (eqs. (|24|) and (|30"])). For the calculation of the CP asymmetry Acp 

Tib — > sgffif) — T(6 — > sqq) . . 

T(b -> s##) + r(6 -> s##) 

the integration over Si and S s should be done. However there are collinear divergences at the 

boundary of the kinematical region. To overcome these divercences we follow the procedure 



given in [14], namely taking a cutoff c in the integration over phase space as: 



^-S s <Si<^(l-c), (19) 



and 

c^<^<^(l-c), (20) 

with c = 0.1. Note that left-right sector gives small contribution to T, however this part is 
responsible for the Aqp- Further Aqp contains only antisymmetric sector. 

3 Discussion 

In the general 3HDM model, there are many free parameters, such as masses of charged and 
neutral Higgs bosons, complex Yukawa matrices, ^- , p^ where i,j are quark flavor indices. 
The additional global 0(2) symmetry in the Higgs flavor space connects the Yukawa matrices 
in the second and third doublet and also keeps the masses of new charged (neutral) Higgs 
particles in the third doublet to be the same as the ones in the second doublet (Appendix 
B). Further, the Yukawa couplings, which are entries of Yukawa matrices, can be restricted 
using the experimental measurements, AF = 2 mixing, the p parameter ]17| and the CLEO 
measurement 



Br{B -> X s7 ) = (3.15 ± 0.35 ± 0.32) KT 4 . (21) 

In our calculations, we neglect all Yukawa couplings except £jy ti ^m,, Pnu an d Pnw by re ~ 
specting these measurements. The same restrictions are done in the model III case and only 
£jy tf and £^ bb are taken into account. 

In this section, we study the s quark energy E s dependency of the differential decay width 
jjjt, differential CP-asymmetry Acp(E s ) and the parameter sinO dependency of the decay 
width T, CP-asymmetry Aqp for the inclusive decay b — > sgg in the framework of the model 
III and 3HDM(02). In our analysis, we restrict the parameters 9, e^ tt and e^ bb (£jv tt and 
£x bb in the model III) using the constraint for \C% |, 0.257 < \C% \ < 0.439 where the upper 
and lower limits were calculated in [IS| following the procedure given in [ PDj] . Here Cy is 



the effective magnetic dipole type Wilson coefficient for b —>■ sj vertex (see |19]|). The above 
restriction allows us to define a constraint region for the parameter e^ tt (£jvtt) i n terms of £® bb 
(^Nbb) an< i ^ i n ^ e 3HDM(02) (the model III). Our numerical calculations based on this 
restriction and the constraint for the angle 9, due to the experimental upper limit of neutron 
electric dipole moment, namely 

<L < 10 _25 e-cm (22) 



which places an upper bound on the couplings with the expression in 3HDM(02) (model III): 
i&^)^ 2 ^ < 1.0 (4® tt S5,)«nfl < 1.0) for m H± « 200 GeV H- Throughout 
these calculations, we take the charged Higgs mass m H ± = 400 GeV, and we use the input 



values given in Table ([!]). 



Parameter 


Value 


m c 


1.4 (GeV) 


m b 


4.8 (GeV) 


A t 


0.04 


m 4 


175 (GeV) 


m w 


80.26 (GeV) 


m z 


91.19 (GeV) 


J^QCD 


0.214 (GeV) 


a s (m z ) 


0.117 


c 


0.1 



Table 1: The values of the input parameters used in the numerical calculations. 



In Fig. 



we plot j^- with respect to the s quark energy E s , for sin6 



and lr t6 | 






< 1. 



dT 
dEs 



is restricted in the region bounded by solid (dashed) lines for Cy > 



0.5,f£ >M = 40m 6 

eff 



(Cj < 0). Dotted line represents the SM contribution. There is a large enhancement in the 
differential decay width for Cj > and C? < in both models, (see [|16| for the model III 
with real Yukawa couplings). In the 3HDM(02), the enhancement is smaller and the restriction 
regions are broader compared to the ones in the model III. 

Fig. | is devoted to the sinO dependence of V for £§ bb = 40 m b and Cj < in the region 
|rt&| < 1. Here F in the model III (3HDM(02j) is restricted in between solid (dashed) lines. 
As shown in the figure, the decay width T increases with increasing sinO. Further, the upper 
and lower bounds for V in the model III exceed the ones in the 3HDM(02) especially for the 
intermediate values of the parameter sinO. Further V can reach the value 10 -3 GeV in both 
models and this is a considerable enhancement compared to the SM one, which is at the order 
of magnitude lO^Gel^ (see ]I3|, JIB]). 

Fig. | (g) shows the E s dependence of A CP (E S ) in the model III (3HDM(0 2 ))- Here solid 
(dashed) lines are the boundaries of the allowed regions of A C p{E s ) for C? > (C? < 0). 
In model III, the restriction region for Acp(E s ) is narrow for Cy > and it has only negative 
values at the order of magnitude 10~ 3 . However, for Cj < 0, this region is broader and 
contains both negative and positive values. The possible values of |A(7p(£' s )| reaches ~ 4% 
for 0.8 GeV < E s < 1.0 GeV. In the 3iLDM (O2), upper and lower boundaries of the allowed 



region for A C p(E s ) are almost coincides for C? > and this region becomes narrower for 
C7 < compared to the one in the model III. In this model \Acp{E s )\ can be ~ 2.5% as a 
maximum value. 

In Fig. |6] and Fig. [7], we represent the sin9 dependence of Acp in the model III and 
3HDM(02)- Acp is restricted in the narrow region bounded by solid lines for Cj > 
and it reaches —0.8 % for sinO = 0.7 in both models. All possible values of Aqp are negative in 
this case. However, for Cj < 0, allowed region becomes broader and Acp can take positive 
and negative values. For this case, |Acp| can reach 3.4%. Note that the restricted regions are 
broader in the model III compared to the ones in 3HDM(02). 

As a conclusion, we get an enhancement in the decay width of the process b —>■ sgg in both 
models, model III and 3HDM(02)- This enhancement is too large to respect the total decay 
width Y tot {b -* sX) = 3.50 ± 1.50 10~ 3 GeV for Cf ff > 0. For C? ff < 0, T can reach the values 
more than two orders of magnitude larger compared to the SM case. Further, we study Acp 
of the process b — > sgg in both models. In the SM, the only source for the CP-violation is the 
complex Cabbibo-Kobayashi-Maskawa matrix elements and almost there is no violation for this 
process. However in the model III and the 3HDM(02), the absolute value of Aqp can reach 
to 3 — 4%, which is a measurable quantity. In addition, we observe that Cj is necessarily 
negative if Aqp has positive values. Therefore the experimental study of the decay width F 
and Aqp of the process b — » sgg can give important clues for the physics beyond the SM and 
also the sign of Cy . 



Appendix 

A The form factors appearing in the b — ► sgg decay 

The function T^ can be divided into two parts: 

rp rpheavy , rplight 

Here T^ avy is the contribution due to the heavy internal quark and neglecting s-quark mass, 
it is given as 

T, v = -i\ t F 2 {( 2p/h ^ah+a^k, _ 2ph ) 

+ — (2cr al3 k°k^g^ + 2 a ua k 2 ^q a -2a fia k 1 t,q a + cr^q 2 ) }m b R . (23) 

where q is the momentum transfer, q = k\ + k 2 , A t is the CKM matrix combination X t = V t bV t * 
and F^hdm j g ^ e form factor 

F 3HDM = p SM {xt) + ^Beyond ^ y[y (24) 

In eq. (p4j) , F 2 (xj) is the magnetic dipole form factor of b —* sg* vertex 
5M -8 + 38x i -39x 2 + 14x; ! -5x; 1 + 18x 2 /nx t 

-^2 ( X <) — ToT i , I \A ' \*o) 



and _F 2 eyon (yt,yi) is the contribution coming from the charged Higgs bosons in 3HDM(0) 

F 2 Beyond (yt, y' t ) = — 2 (C N u a + £vi§) (&,« + &,Jr) G i (*) 

m t v ts V tb 

1 - V* - V 

+ (Cn,U + CnJcTT^) i^N,bb + £,N,sb71~)G2(yt) 

m t m b Vts ' v tb 



2- 



~~2 \PN,tt + Pv,ic T a* ) \PN,tt + PN,tc-ir , 

™t v ts v tb 



with 



+ —(ffiu + ffitM)(pN& + f%Jr)G2(3/ t ) (26) 



Gi(y) = u J 1+ ((-l + y)(-2-5y + y 2 )+6ylny)), 

G 2 (y) = 2( J +y)A (y(3-4y + y 2 ) + 2(-l + y)ylny). (27) 



where yt = m * and y' t = -^- (see appendix B). In eq. ( p6[) we used the redefinition 



HP) U '° = ^ apf' D ■ (28) 

In eq. (|To|) , T l ^ ht contains two different parts related to the light internal quarks. The 
first one is obtained from T^ avy with the replacement F 2 3HDAI _^ _pSM^ an( ^ ^e second one 
T 2 l ^ v is the contribution due to the non-vanishing kl xternal /mf ight terms 

- 2AF 1 (j v k 2li -^k lv )L} (29) 



where 



with 



9 3 z 3 

A- 5 9 0-(*) , 8Q (s) 2 

9 Z 6 Z 6 

Ai 5 = -1-2^^, (30) 



Qo(-z) = -2 - (w+ - w_)(£n hi7r) , 

Q_(z) = -(Zn— + m) 2 . (31) 

2 Mi 



Here the parameter w± is 



and 



u ± = ^(l±Jl-- z ), (32) 



g 2 
z = — =■, i = u,c . (33) 

Finally the function T^ reads as 

T M „ = a R (T^ v + T h ^ y {F! HDM - -if M (0))) + a L ig? . (34) 

Here or, «£ are real parameters to seperate the parts with factors R and L in the fuction T^. 
With this parametrization T can be written as 

T = a R T R + a\Y L + a R <x L T LR \ ah ^ aR ^ x . (35) 

Note that the expressions for the model III case can be obtained by disregarding the Yukawa 

{D. 
ij 

10 



i- -U,{D) . 

couplings p N \ • in eq 



B 3HDM(0 2 ) 



In the multi-Higgs doublet (n > 2) models, the Higgs sector is extended and therefore the 
number of free parameters, namely, masses of charged and neutral Higgs particles, Yukawa 
couplings, extremely increases. In our problem we choose n = 3 and to overcome the difficulty 
coming from the large number of free parameters we consider 3 Higgs scalars as orthogonal 
vectors in the Higgs flavor space, denoting by the index to", where m = 1,2,3. At this stage 



we introduce a new global 0(2) symmetry on the Higgs sector [22 



0i = 01 , 

2 = cos a 4>2 + sin a 03 , 

3 = —sin a 02 + cos a 03 , (36) 

where a is the global parameter, which represents a rotation of the vectors 02 and 03 along 
the axis that 0i lies, in the Higgs flavor space. This transformation keeps the kinetic term of 
3HDM Lagrangian invariant: 



C-Kinetic — {D^4>i) D 



V, 



where 



(&*& - i 9 -B^ - i\4*\w») (37) 



1,2,3. (3* 



12 



The invariance of the potential term 

V{<f> x , 02,0 3 ) = c 1 (0+0 1 -^/2) 2 + c 2 (0+02) 2 

+ 03(0^03)' + C 4 [(0^0l - V 2 /2) + 0+02 + + <p,| 

+ C 5 [(0 + 1 )(0 + 02)-(0^02)(0 2 f 0l)] 

+ C 6 [(0 + 1 )(0+0 3 )-(0+03)(0 3 + 0l)] 

+ C 7 [(0+0 2 )(03 f 03)-(02 f 03)(03 + 02)] 

+ c 8 [i?e(0+0 2 )] 2 + c 9 [i?e(0+0 3 )] 2 + c lo [ J Re(0 + 3 )] i 

+ c 11 [/m(0 + 02)] 2 + c 1 2[/TO(0 + 03)] 2 + c 13 [/m(0 + 03)] 2 + c 1 4 (39) 



11 



12 



can be obtained if the following conditions on the free parameters are satisfied: 

c 5 = c 6 ) c 8 = c 9 > c ll = c 12 j 

c 2 = c 3 = c 7 = cio = . (40) 

This implies that the masses of new particles are the same as the older ones, namely, 

v 2 
m F± = m H± = c 5 — , 

v 2 
m H 3 = m H i = c 8 — , 

,2 



V 

m H 4 = m H 2 = en— , 



(41) 



Further, the application of this transformation to the Yukawa Lagrangian (eq.([|)) keeps it 
invariant if the transformed Yukawa matrices satisfy the expressions 

pU{D) ±U(D) JU(D) ■ 

iij = iij cos a + p^ 'sin a, 

-IU{D) FU(D) ■ -(7(D) , An \ 

p y v = -&/ 'sma + Pif 'cos a. (42) 

and therefore 

/g/U{D)^+g/U(D) + ^U(D)^+-,U(D) _ ^U(D)j+£U(D) + (-U(D)y-U(D) ^ /^ 

which allows us the following possible parametrization of the Yukawa matrices ^ U( - D ^ and p u ^: 

£U(D) =e U(D) cose ^ 

p u = e u sin 9 , 
p D = ie D sin9, (44) 

where e u ^ D ' are real matrices satisfy the equation 



Here T denotes transpose operation. In eq. (|44]) , we take p D complex to carry all CP violating 
effects in the third Higgs scalar. 

Therefore we can reduce the number of free parameters taking the new charged and neutral 
boson masses as the same as the older ones and connecting the Yukawa matrices ^^ and 
pU{D) us j n g thg expression eq. (|45f) . 

12 



Note that, neglecting the off-dioganal Yukawa couplings, the expression for F 2 eyon (y t , y' t ) 
(eq. (p6|)) can be written as 

F Beyon d{yt) = J_ {f f Gx{yt) + J_f g£ >bb G 2 (y t )(cos 2 9 + isin 2 9) (46) 
mi rntrrib 



13 
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Figure 1: -4Sr as a function of E s for fixed £nm> = 40 m;,, sm# = 0.5 and \r t b\ = \jir^\ < 1 in 
the model III. Here jjjr is restricted in the region bounded by solid (dashed) lines for Cj > 
(Cj < 0). Dotted line represents the SM contribution. 
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Figure 2: The same as Fig. | but for 3HDM(0 2 ). 



Si 
(-1 




Figure 3: T as a function of sin # for C? < 0, £§ bb = 40 m&, and \r tb \ < 1. Here T is restricted 
in the region bounded by solid (dashed) lines for the model III (3HDM(02))- 
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Figure 4: Acp(E s ) as a function of E s for fixed £^ feb = 40 m^, sinfl = 0.5 and |r t &| < 1 in the 
model III. Here A C p(E s ) is restricted in the region bounded by solid (dashed) lines for C? > 
(C 7 e// < 0). 




Figure 5: The same as Fig. |, but for 3HDM(0 2 
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Figure 6: Aqp as a function of sin 8 for £ 
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N,bb 



AOrrib and Ir^l < 1 in the model III. Here Aqp 



is restricted in the region bounded by solid (dashed) lines for C? > (C- < 0) 
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Figure 7: The same as Fig |6], but for 3HDM(C>2). 
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